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Abstract 

We present a detailed numerical calculation of the upper critical field H C 2(T) 
for a bulk extreme type-II superconductor. Particular emphasis is placed on 
the high-field, low-temperature regime of the HT-phase diagram. In this 
regime it is necessary to go beyond the standard semi-classical theory and 
include the effects of Landau quantization of the electronic motion on the 
superconducting state. The presence of Landau level quantization induces an 
upward curvature in H C 2(T) at ~ 10% of T c q for those superconducting sys- 
tems in which the slope of H&iT) at T c q is > 0.2 Tesla/Kelvin. We construct 
a simple analytical model that can account for this behavior based on the 
renormalization of the BCS coupling constant by the off-diagonal pairing of 
electrons on Landau levels. 
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There has been considerable recent evidence of an anomalous "divergence" of the up- 
per critical magnetic field H c2 (T) at low temperatures in a number of "low-temperature" 
high- temperature superconductors (HTS). Upward curvature and the anomalous divergence 
was observed in B^S^CuOj, [|TJ and T^E^CuOq thin films 0; K z Bai_ x Bi03 single crys- 
tals ; Tl 2 Mo 6 Se 6 compound |5[] as well in borocarbide intermetallic superconductors 
||. These observations are in contradiction with the standard semi-classical Werthamer- 
Helfand-Hohenberg (WHH) theory |7]||] which predicts the saturation of H C 2(T) at low tem- 
peratures. WHH theory yields a finite value of the upper critical field at zero temperature, 
H^ HH (0), and offers an elegant method for determining this field in clean superconductors 
as 0.693T c0 (— dH c2 /dT) T=Tc0 , where (dH/dT) T=Tc0 is the slope of the upper critical field eval- 
uated at the zero-field transition temperature T c0 (later in the text called the WHH-slope). 
For almost three decades, the WHH method has enabled experimentalists to determine 
H C 2(T) at low temperatures from the high temperature and low magnetic field data (around 
T c0 ) in conventional low-temperature superconductors . Recent advances in obtaining high 
magnetic fields in a laboratory environment || have revealed deviations of the experimental 
data from WHH theory at low temperatures. 

WHH theory is based on the semiclassical phase-integral approximation originally due 
to Gor'kov [|K|] which depends on the following assumption: the bending of the semiclassical 
paths of the electrons by the magnetic field is negligible over the range of the single-particle 
Green's function at zero magnetic field. The latter is given by Vp/i^l-irksT), while the radius 
of a semiclassical path equals l 2 kp, where vp and kp are the Fermi velocity and wave vector, 
and I = ^Jhc/eH is the magnetic length. For a clean superconductor, this assumption reads 
l 2 hp ^> Vp/finksT) or equivalently huj c <C ksT, where uo c = eH/m*c is the cyclotron 
frequency. In dirty superconductors (with significant impurity concentration), it translates 
to Huj c <C T where T is the scattering rate due to disorder. From these conditions, it 
is clear that the semiclassical WHH theory ignores quantum Landau level (LL) effects, an 
assumption justified at low fields where the electrons occupy a huge number of closely spaced 
Landau levels (LL's). In this case, the temperature and/or impurity scattering broaden LL's 
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and reduce the significance of LL quantization. In the opposite limit of high magnetic fields 
and low temperatures, where k^T <C hu c and T -C hu c , the quantization effects are of crucial 



importance. It was shown by Tesanovic et al. [TT| that the discreteness of the LL's leads to 



a breakdown of the semi-classical picture. When the LL level structure is fully accounted 
for in the BCS theory, one finds that, for a pure case or for a moderate level of impurities, 
superconductivity does not terminate above the semiclassical H^ HH (T) line, but exhibits 
reentrant behavior where superconductivity is enhanced by the magnetic field H (where 
H 3> H^ HH (T)) fll,ll|. A similar reentrant phase at high fields was recently proposed 



in organic quasi-one-dimensional superconductors JT^]. As a consequence of the underlying 
LL structure H c2 (T), or rather T C (H), develops oscillations near H^ HH (0) signalling the 
passage of the LL through the chemical potential. Similar types of quantum oscillations, 
with the same origin, have been predicted in various other measurable quantities which are 



particularly pronounced in two-dimensional systems JT2|.[T4 . 

The LL quantization of electronic orbits within the superconducting state is now an ex- 
perimental fact in numerous extreme type-II systems. The strongest evidence so far comes 
from recent observations of de Haas-van Alphen (dHvA) oscillations deep in the mixed 
state of various superconducting materials ranging from A-15's and boro-carbides to high-T c 
cuprates [p!q-|17| . What all these materials (some of them also being samples with a "diverg- 
ing" upper critical field) have in common is their extreme type-II character: the WHH slope 
in these systems is comparable to ~ 0.2 Tesla/Kelvin and is often higher In such sys- 

tems the cyclotron splitting of LL's near H^ HH (0), hu c2 (0), where o; c2 (0) = eH c2 (0)/m*c, is 
comparable to ksT cQ and there is a large region in the H-T phase diagram in which the LL 
structure within the superconducting phase is well defined whith hu c > A(T, H), ksT and T 
(here A(T, H) is the BCS gap) |Tj|. The boundaries of this high field and low-temperature 
region in the H — T diagram, H* and T*, extend to fields as low as H* ~ 0.5H C 2(0) and 



temperatures as high as T* ~ 0.3T c o ||16|| . In contrast, the size of this region in conventional 
type-II superconductors (like Nb) is negligible. Within the BCS theory, the scale of the 
cyclotron splitting between LL's near H C 2(0) in conventional systems is set by the condensa- 



tion energy, ~ (ksTco) 2 /Ep, and should be much smaller than either the thermal smearing, 
~ fcfiT, or the BCS gap A(T,H). Additional smearing due to the disorder T makes this 
high-field and low-temperature region in H-T diagram irrelevant. 

In this paper we present a detailed numerical calculation of the upper critical field H C 2(T) 
for the extreme bulk type-II superconductor with particular emphasis on the high field and 
low temperature region in the H-T diagram. We assume a simple isotropic model based on 
the mean-field theory on LL's developed in our previous work fl6|]. We find LL structure 
inducing upward curvature in H C 2(T) at ~ 10% of T c q in these systems. In order to account 
for such behavior, we construct a scheme based on the renormalization of the BCS coupling 
constant by the off-diagonal pairing of the electrons on LL's. The model presented in this 
paper can be used to describe behavior of a cubic material such as K x Bai_ :r Bi03 0,f|. 
Furthermore, it can be extended to those anisotropic systems in which the hopping energy 
t 1 - is larger than a cyclotron gap at high fileds, i. e. t > hu c . In such case, an isotropic 
model, with renormalized masses, correctly captures high-field behavior of a system. In the 
opposite limit of a highly anisotropic material or a quasi-two-dimensional system, mean-field 
theory is of a limited value due to strong fluctuations [ |TTf . 

Within mean-field theory, the transition line T C (H), for any field H, is described as the 
solution of the self-consistent equation 

A(r) = ^£F(r,r;^) (1) 

when the amplitude of the superconducting order parameter A(r) goes to zero. V is the 
usual BCS interaction strength and (3 = 1/ksT. F(r, r; uo) is the anomalous Green's function 
with uj = 2irkBT(m+ 1/2) being the electron Matzubara frequency (the Matzubara index m 
is dropped in order not to be confused with the LL index later in the text). This anomalous 
Green's function is constructed from the wavefunction of the individual electrons in a mag- 
netic field (for review see reference |19[]). For a clean three-dimensional (3D) superconductor 
in high magnetic field equation (ffl) reduces to the form 
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V ~ (n + n')\ fdk z ^ 1 1 ( . 



2nP(3 n f^ =Q n\n'\2 n + n ' J 8tt ^ tcu - e n± (k z ) -ilo - e n ^(k z ) 

where n and n! are indices of the Landau levels participating in the superconducting pairing 
and k z is the momentum along the field direction. The electronic energies in the magnetic 
field 

tn±{k z ) = H + hu c (n +\)T9^-, = ^(kl- k Fn± ) (3) 

are measured from the chemical potential \i. In the presence of Zeeman splitting (5 ^ 0), 
for each value of n and for each index "plus" and "minus" there is a corresponding Fermi 
momentum kp n ± determined by the conditions (for T///<1) 



2m* HeH 
k Fn ± = )l-^ r (ii-huj c (n + 1/2) ±g— ) (4) 

n e=^pH k Fn± (5) 

where n e is the electronic density and the sum is over all real kp n ±. In deriving self-consistent 
equation |2| we assume that the order parameter A(r) has a standard Abrikosov form [30 



In the LL representation Abrikosov order parameter belongs to the lowest Landau level 
(LLL) for the Cooper charge e* = 2e having the lowest kinetic energy of the center-of-mass 
motion of Cooper pairs and therefore the highest transition temperature. Note though, that 
in two dimensions (2D) special circumstances might arise where the contributions to the 
order parameter from the higher LL's become competitive leading to the higher transition 
temperature However, in 3D, higher LL contributions lead to the significantly lower 



transition temperature 11,21 



In the framework of the BCS theory, the density of states (DOS) of the electrons in 
a magnetic field exhibits a divergence whenever the bottom of a LL crosses the chemical 
potential. As a consequence, strong oscillations in T C (H) develop at high fields and low tem- 
perature in a clean superconductor |TT[. These oscillations, in principle, could be observed 
in extremely clean samples which is rarely the case in an experimental setup. Therefore, we 
consider a more realistic system for which the presence of the impurities and imperfection 



(or disorder in general) is taken into account. We assume that the disorder in the sample 
leads to isotropic broadening of the LL's, the size of which is measured by T = h/2r , where 
l/2r is the scattering rate due to the disorder. As long as Y/ftio c <C 1 the discreteness of 
the LL structure is preserved but the quantum oscillations are greatly reduced. In this way 
we are left only with the task of numerically examining the overall rising trend of T C (H) at 
high fields. The broadening of the LL can be most easily included into the self-consistent 
equation @ by a substitution iu — > iu — T in the Matzubara frequencies. 

After integration over the momenta k z and summation over the Matzubara frequencies 
(with the Debye frequency Q as a UV cut-off), the equation (0) can be put in the form 

1 - (, + ,-)! fajjl r _ _ 



A n y =0 n\n>\2 n + n ' 2tt1 2 N 3D {0) 4 



f2 ,hu c (n — n' ± g/2) 

2^1+ = h z 

* 2vrT 4ttT 



1 .hw e (n-n'±g/2) 
Zi± = - 2 +* ^ (6) 

where for each value of the index n there are two terms, one for the "+" and one for the "-" 
sign. A = KiV3£)(0) is the BCS coupling constant with N 3 d(0) being the single-spin electronic 
DOS in zero field calculated at the Fermi energy Ep. *&(z) is the Digamma function of the 
complex variable z, and z* denotes complex conjugation. Nicm(0) /2nl 2 is the DOS of the 
center-of-mass of the Cooper charge at the chemical potential and it compares to N^d(0) as 

N 1CM {o) r 



2nl 2 N 3D (0) 2 5 / 2 t/Ef7kJ~ c (C 2 + (^) 2 ) 1/2 ( X /C 2 + (-^) 2 - Q 1 / 2 



(7) 



where 



= jj_ _ n + n' + 1 

huj c 2 K ' 



Our main goal is to numerically solve the self-consistent equation (§) for the transition 
temperature T C (H) taking realistic values for the coupling constant A as well as moderate val- 
ues for the disorder parameter T /Ep . In particular, we want to incorporate in our calculation 
the crucial observation that the WHH-slope (—dH c2 / (1T)t=t c0 reported in the experiments 
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|T|-|| is always much larger than pa 0.2 Tesla/K. Therefore, we chose a coupling constant 
A so that the WHH-slope when expressed in the dimensionless units [— du c / d(k B T)] T=Tc0 is 
larger than 0.27 0. 

Figure 1 shows the numerical solution of equation (|6|), i.e. the plot of the upper critical 
field H c2 (T) rescaled by H^ HH (0) vs. T c /T c0 for the model system with coupling constant 
A = 0.35 (full line) and A = 0.4 (dotted line), the disorder parameter T/Ep = 0.025 and zero 
Zeeman splitting (g = 0). For the model with A = 0.35 the WHH-slope in dimensionless units 
is [— dui c / d{kBT)]T=T c0 = 0.8 which corresponds to (—dH C 2/dT)T=T c0 — 0.58 Tesla/K. When 
A = 0.4, the WHH slope in dimensionless units is 1.108, or 0.813 Tesla/K in conventional 
units. The upper critical field starts to deviate from the one predicted by the WHH theory 
at T c /T c0 pa 0.1 exhibiting an anomalous divergence at low temperature as predicted by 
Tesanovic et al. [|Il|] and recently seen in experiments While disorder completely 

washes away the quantum oscillations of the upper critical field at low temperatures for 
A = 0.35 case, they are still visible for the model with the A = 0.4. 

In Figure 2 we plot the ratio huo c / Ep (which is proportional to H c2 , the proportionality 
constant being set by the material properties of the particular superconductor) vs. T c /T c0 for 
the model system with A = 0.35 in order to illustrate the effect of Zeeman splitting (g ^ 0) 
on the upper critical field. The overall tendency of the Zeeman splitting is to suppress T c in 
the region of the H-T diagram where /cgT < hu c , while the transition line is not affected by 
Pauli pair breaking in the high temperature-low field portion of the diagram. An interesting 
situation arises when effective g factors are very close to g = 2m/m*, g = 4m/m* or any 
other even integer (we take the effective cyclotron mass m* of the order of the electron mass 
m further in the text). Unlike the odd g factor case (g = 1 is an example shown in Figure 
2), the anomalous divergence of the upper critical field at low temperature is not destroyed 
by Zeeman splitting for even integer g factors. The temperature at which H c2 (T) starts 
to deviate from the WHH line is reduced by only few percent from the g = case. This 
situation can be understood as follows: When g = 2 or g = 4 the Zeeman splitting is equal 
to the cyclotron splitting making the nth spin-up LL degenerate with the (n + l)th (for 



g = 2) or (n + 2)th (for g = 4) spin-down level. In these cases the off-diagonal terms in 
equation (||), describing LL's separated by the cyclotron gaps, become effectively diagonal 
(i.e. degenerate) leading to the anomalous divergence in H c2 (T) (see the discussion below). 
When g is an odd integer or a fraction, the LL spin degeneracy is completely lifted resulting 
in the suppressed H c2 (T) at low temperatures with a downward curvature. 

The upward curvature of the upper critical field H C 2(T) at temperatures T < 0.1T c o is 
the consequence of the competing tendencies of diagonal, n = n', and off-diagonal, n ^ n', 
terms in the self-consistent equation @. The diagonal terms correspond to the Cooper pairs 
formed by electrons in the same LL while the off-diagonal terms represent the electronic 
pairing of the LL's separated by hu c or more. Only the diagonal terms posses Cooper 
singularity and therefore lead to the increasing trend of the transition temperature as a 



function of field pl|J12|Jl^] . In clean systems, they also produce strong quantum oscillations 
in T C (H) for fields H > H^ HH (0). Ultimately, the diagonal terms in (|6|) lead to the reentrant 
superconductivity for magnetic fields H 3> H^ HH (0). This Quantum Limit regime, reached 
when all electrons in the system occupy a single, lowest Landau level, was predicted and 



investigated by Tesanovic et al. [IT] but is not of interest in our calculation since it can 
be realized only in extremely high fields. The off-diagonal terms in (|6[) start to dominate 
at lower fields H m H^ HH (0) creating a counter-effect to the rapidly decreasing diagonal 
terms. These non-singular terms lead to the smooth crossover from the diverging high-field 
transition line to the low-field WHH line. This crossover results necessarily in the upward 
curvature of the upper critical field H c2 (T) at low temperatures. Our goal in this paper is 
to closely examine this crossover behavior of the off-diagonal terms and to account for the 
resulting upward curvature in a simple analytic model. First, we notice that the effective 
role of the off-diagonal terms in the self-consistent equation (|6]) is to renormalize the BCS 
coupling constant A into a new, field and temperature dependent, constant X(H,T) through 
the substitution 1/A -> 1/A - 1/\'(H,T) = 1/\{H,T) where 1/\'(H,T) accounts for the 
off-diagonal, n ^ n', terms . The size of the off-diagonal terms in (|6|) grows as the magnetic 
field is lowered, leading to the effective increase in the renormalized BCS coupling constant 
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X(H,T) rendering a non-zero transition temperature at lower fields. 

With this approach the solution of the self-consistent equation can be written in analyt- 
ical form as 

-i" 



T C {H) = 1.134ft exp 



2irP 



^ N ln (0)(2n)\ 



(9) 



X(H,T C ) 

where iV\ n (0) is the ID density of states given by (|?p and (|^), where n = n' . When 
2j[ksT < hu c the contribution of the off-diagonal terms 1/X'(H,T C ) to the renormalized 
coupling constant can be obtained by systematic expansion in 2irkBT c /huj c . Assuming that 
the number of the occupied LL's n c = Ep/hu! c is much larger than one, i.e. that the sys- 
tem in question is far away from the Quantum Limit regime (in a typical experimental 
setup this is indeed the case, since usually n c > 30), and ignoring the quantum oscillations 
(they are damped by disorder), we perform the expansion of the off-diagonal terms up to 
(2iikbT c /hu! c ) 2 . Then, it is possible to rewrite the equation (||) in very simple analytic form 



as 



T C (H) = 1.134ft exp 



^ \X X'(H,T C ) 



(10) 



where 



1 



X'(H,T C 



( 

exp + exp 



71 r 



+ 



2^Tvn c 



if 2i\kBT c 
In Jn c H 



2n c \ huj c 



2 1 

+ 



7T~ 



— In (n c ) 



' 2nk B T c 



(11) 



Equation ([!(]) can be solved by iteration or can be used as a formula to fit the experimental 
data. Figure 3 shows T C (H) computed from (|9p with 1/ X'(H,T C ) in ([TT]) expanded to leading 
order (dotted line) and to second order (dashed line) in 27ikBT c /fuj c compared to the exact 
numerical solution of self-consistent equation (^j for A = 0.35. The agreement is excellent 
in the region where 2irkBT < hu c , where the expansion ( |11~D is valid. The simple analytic 
form ( [TDD accounts very well for the diverging upper critical field H c2 (T) and its anomalous 
upward curvature. When 2ixksT = Tioj c this expansion breaks down as indicated in Figure 
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3 by the straight line. At low fields and high temperature, where 27rfceT > hu c , there is a 
large deviation of T C (H) from the exact numerical solution of signalling the breakdown 
of the expansion. 

To conclude we have presented a detailed numerical calculation of the upper critical field 
H c2 (T) for a three dimensional extreme type-II superconductor characterized by a large 
WHH slope. We find that the Landau level quantization induces an upward curvature in 
H C 2(T) at temperatures ~ 0.1T c o. We account for this behavior through renormalization 
of the BCS coupling constant A by the off-diagonal pairing of the electrons in the Landau 
levels. Our work, based on the simple BCS model, reproduces qualitatively observations of 
the anomalous behavior of the upper critical field in the experiments ||^[], but it cannot 
account quantitatively for the large deviations in H c2 (T) from the WHH line for temperatures 
higher than predicted by our theory. We believe that extension of this work to a more realistic 
strong-coupling model within the Landau level framework will improve agreement with the 
experimental results. 

We are grateful to Prof. Zlatko Tesanovic for many useful discussions. This work has 
been supported by NSF Research Opportunity Award DMR-9415549, Amendment No. 003. 
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FIGURES 

FIG. 1. Plot of the upper critical field H C 2 computed from (6) vs. T c /T c q for a 3D supercon- 
ductor with the BCS coupling constant A = 0.35 and A = 0.40 and with no Zeeman splitting. H C 2 
is rescaled by H% HH (0) = 0m3T c0 (-dH c2 /dT) T=TcO . 

FIG. 2. Upper critical field H C 2 ~ Tilo c /Ef for a 3D superconductor with non-zero Zeeman 
splitting. 

FIG. 3. Comparison of the exact numerical solution of (6) (full line) and the solution of (9) 
obtained by the expansion of the off-diagonal terms in (6) up to the leading (dotted line) and the 
second order (dashed line) in 27r&BT c /ftu; c ). To the left of the straight line 27t/cbT c = huj c this 
expansion is valid. 
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